DUAL LUKACS REGRESSIONS FOR NON-COMMUTATIVE VARIABLES 



KAMIL SZPOJANKOWSKI AND JACEK WESOLOWSKI 

Abstract. Dual Lukacs type characterizations of random variables in free probability are studied here. 
First, we develop a freeness property satisfied by Lukacs type transformations of free-Poisson and 
free-Binomial non-commutative variables which are free. Second, we give a characterization of non- 
commutative free-Poisson and free-Binomial variables by properties of first two conditional moments, 
which mimic Lukacs type assumptions known from classical probability. More precisely, our result is a 
non-commutative version of the following result known in classical probability: if U, V are independent 
real random variables, such that E{V{1 - U)\UV) and E{V^{1 - Uf\UV) are non-random then V has 
a gamma distribution and U has a beta distribution. 



1. Introduction 

Characterizations of non-commutative variables and their distributions is a field which develops through 
non-commutative probability with results which parallel their classical counterparts. It is not completely 
well understood why the results mirror so much these from the classical setting since the nature of objects 
under study seems to be much different. 

An example of such a result is the Bernstein chracterization of the normal law of independent random 
variables X and Y by independence oi X + Y and X — Y in classical probability |3j (see also [21]), and a 
characterization of non-commutative semicircular variables X and Y which are free and such that X+Y and 
X — Y are free by Nica in [^H] ■ Similarly, the classical characterization of the normal law by independence 
of the mean X = ^ J27=i -^i^ ^^'i empirical variance = S"=i i-^i ~ : where (Xi)i=i_..._„ are 
independent, identically distributed real random variables from [35] is paralleled by a non-comutative 
charactrization of the Wigner law exploiting freeness of X = i Y^^=i ^^'^ — TT^T ^"=1 O^i ~ ■^)^ 
built on free indentically distributed non-comutative random variables {Xi)i=i,...,n - see [20] . 

In this paper we are concerned with the celebrated Lukacs characterization of the gamma distribution, 
[2 7) . It says that if X and Y are positive, non-degenerate and independent random variables and such 
that 

(1) U = — ^ and V = X + Y 

are independent then X and Y have gamma distributions, G{p,a) and G{q,a). Here by the gamma 
distribution G(r, c), r, c > 0, we understand the probability distribution with density 

f{x) = — a;'-le-^"/(o,oo)(a;)• 
The direct result: li X ^ G{p, a) and Y ^ G{q, a) are independent then U and V, defined through ([T]), 
are independent; is rather simple. It suffices just to compute the jacobian of the bijective transformation 
(0, oo)^ 9 {x,y) e (0,1) X (0,oo) and to follow how the densities transform. Immediately it 
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follows also that V ~ G{p + q, a) and [/ is a beta random variable Bi{p, q), which has the density 

The same computation while read backward proves the opposite implication: if U and V are independent, 
U ^ G{p + q, a) and V ^ Bj{p, q) then X = UV and Y = {1 — U)V are independent, X ^ G{p, a) and 
r ^G(q,a). 

For random matrices the role of the gamma law was taken over by Wishart distribution defined, 
e.g. on the cone V_|_ of non-negative definite real n x n symmetric matrices by the Laplace transform 

-^(^) ~ ( dct°a+s) ) positive definite a and p € {0, 5, |, . . . , U (^^, 00), and for s such that 

a + s is positive definite. If p > then Wishart distribution has density with respect to the Lebesgue 

measure on V+ of the form 

/(x) cx (det x)P-'^ e" '^''/v+ (x). 

Matrix variate beta distribution, in the case of real n x n matrices, is a probability distribution on the 
set 2? = {x e V+ : I — x e V+} defined by the density 

(7(x) cx (det x)P-i (det (I - x))«-\ 

where the parameters p,q G {^^^, 00). 

Analogues of Lukacs characterizations have been studied since Olkin and Rubin [211 [32] • In particular, 
Casalis and Letac [TS] obtained such a characterization in a general setting of probability measures 
on symmetric cones (including positive definite real symmetric and hermitian matrices), though they 
assumed an additional structural invariance property. Under smoothness conditions on densities Bobecka 
and Wesolowski [4] proved that if X and Y are independent V+-valued random matrices and U = 
(X + Y)~ 2 X (X + Y)~2 and V = X + Y are independent then X and Y are Wishart matrices. For 
recent extensions see [H |8] . 

In the context of Lukacs type characterizations of distributions of random variables in non-commutative 
setting the analogue of gamma distribution is free-Poisson (Marchenko-Pastur) distribution (Note that 
this analogy follows neither Berkovici-Pata bijection [2] , nor analogy between classical and free Meixner 
distributions defined in jlj and developed in [9] i.e. it is not free gammma distribution). Let X and Y 
be free non-commutative variables having free-Poisson distributions (with properly defined parameters). 
Define 

U = (X + Y)-i/2x(x + Y)-i/2 V = X + Y. 

One would suspect that by the analogy to the classical case or to the matrix variate situation, U and 
V are free. This is still an open problem. The closest result has been derived in [14] (referred to by CC 
in the sequel). They proved that for complex Wishart independent matrices X and Y the matrices U 
and V defined as for the real case above are asymptotically free and the limiting (in non-commutative 
sense) distributions of U and V were derived to be free-Poisson and a distribution which, by the analogy 
to classical (univariate or matrix- variate) cases could be called free-beta, but it has already been known 
under the name free-binomial (for details, see Sect. 7 in CC; consult also the first part of Section 2 below 
where a complete description of the set of parameters of this distribution is presented) . The main result 
of Section 3 goes in the opposite direction. Assuming that free variables U and V have, respectively, 
free-binomial and free-Poisson distributions we prove that X and Y are free with suitable free-Poisson 
distribution each. This is done through developing some ideas from CC. 

The direct non-commutative version of Lukacs characterization, saying that if X and Y are free and U 
and V, as defined above, are free then X and Y are free-Poisson was obtained in Prop. 3.5. 

The classical Lukacs characterization can be obtained with weaker assumptions than independence of 
U and V. Such assumptions may be formulated in the language of constancy of regressions. For instance, 
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it is known that if X and Y arc positive, non-degenerate and independent and 

(2) ¥.{X\X + Y) = c{X + Y) and ¥.{X^\X + Y) = d{X + Yf 

for some real numbers c and d then X and Y are necessarily gamma distributed, G{p,a) and G{q,a), 
where the parameters p and q depend on c and d. This can be traced back to Bolger and Harkness, 
[S]. But the result is also hidden as one of special cases in the celebrated Laha and Lukacs paper 
Regression versions of Lukacs type characterizations of Wishart random matrices were obtained in |24) 
and more recently in the framework of regressions of quadratic forms in [25j [26] . 

The non-commutative version can be found in [9], which is devoted mostly to Laha-Lukacs type 
characterizations. They assumed that (/3(X|X -I- Y) = i(X -I- Y) which is an analogue of first part of ^ 
for identically distributed X and Y, but instead of a direct non-commutative version of the second part 
of ^ they considered, following the classical setting of [23] a more general condition 

ip{X^\X Y) = a(X + Yf + b(X + Y) + cl. 

They used a cumulant approach to derive possible distributions of X and Y. A related results in converse 
direction for free-Poisson variables has been given recently in [TH] 

Our aim here is to consider the dual regression scheme of Lukacs type. In the classical setting it means 
that we take idependent U which is (0, l)-valued and V which is positive and assume that 

E{{1 ^ U)V\UV) ^ c and E{{1 - U)^V^\UV) ^ d 

for some constants c and d. It was proved in ^ that then necessarily U and V are, respectively, beta 
and gamma random variables (see also [11] for a more general characterization). A version of this 
characterization in the non-cummutative setting which is considered in Section 4, is the main result of 
this paper. 

Next section is devoted to introduce basics of non-commutative probability we need to explain the 
results and derivations. 

2. Preliminaries 

Following [37j or |30j we will recall basic notions of non-commutative probability which are necessary 
for this paper. 

A non-commutative probability space is a pair (A, ip), where A is a unital algebra over C and ip : A ^ C 
is a linear functional satisfying (f){T) — 1. Any element X of ^ is called a (non-commutative) random 
variable. 

Let H he a Hilbert space. By B{H) denote the space of bounded linear operators on H . For A C B{H) 
and tp E H we say that {A, ip) is a VF*-probability space when A is a, von Neumann algebra and ip 
is a normalized, faithfuU and tracial state, that is \\(p\\ = 1, <(5(X^) {'K'^p,tp) = iff X = and 
ipiXY) = ip{Y* X) for any X,Y e B{H). 

The ^-distribution /i of a self-adjoint element X e ^ C B{H) is a probabilistic measure on M such that 

p{W) ^ [ t'^i{dt) Vr 1,2,... 

In a setting of a general non-commutative probability space {A, p), we say that the distribution of the 
family 0^i)i=i^...^q is a linear functional ^Xi....,x, on the algebra C {xi, . . . ,Xg) of polynomials of non- 
commuting variables Xi, . . . , Xq, defined by 

MXi,...,x,(P) = ^{Pi^u ■ ■ • ,X,)) VP € C . . . 

Unital subalgebras Ai C A, i ~ l, . . . ,n, are said to be freely independent if </?(Xi, . . . ,Xfc) = for 
Xj G Ai(^j-) , where € {1, . . . , ti}, such that ipi^^j) = 0, j = 1, . . . , A;, if neighbouring elements are from 
different subalgebras, that is z(l) ^ i{2) ^ . . . ^ i{k—l) ^ i{k). Similarly, random variables X, Y e ^ are 
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free (freely independent) when subalgebras generated by (X, I) and (Y, I) are freely independent (here I 
denotes identity operator). 

For free random variables X and Y having distributions /i and v, respectively, the distribution of X + Y, 
denoted by ^ ffl j^, is called free convolution of /i and v. 

For self-adjoint and free X, Y with distributions /i and ^, respectively, and X positive, that is the 
support of /I is a subset of (0, cxj), free multiplicative convolution of /i and v is defined as the distribution 
of VX YVX and denoted by ^Mv. Due to the tracial property of ip the moments of Y X, X Y and \/X Y\/X 
match. 

Let X = {Bij i?2, • • •} be a partition of the set of numbers {1, . . . , fc}. A partition x is a crossing 
partition if there exist distinct blocks € x ^^id numbers ii,i2 G Br, ji,j2 € Bs such that 

ii < ji < 12 < j2- Otherwise x is called a non-crossing partition. The set of all non-crossing partitions 
of {1, . . . , fc} is denoted by NC{k). 

For any fc = 1,2,..., (joint) cumulants of order k of non-commutative random variables Xi, . . . , X„ 
are defined recursively as fc-linear maps TZ^ : C (a;^, j = 1, . . . , fc) — >■ C through equations 

(^(Yi,...,Y„,) = W ^^\B\{x^, I e B) 

xeNCim) Bex 

holding for any Y.; S {Xi, . . . , X„}, i = 1, . . . , m, and any to = 1, 2, . . ., with \B\ denoting the size of the 
block B. 

Freeness can be characterized in terms of behaviour of cumulants in the following way: Consider unital 
subalgebras {Ai)ii£i of an algebra ^ in a non-commutative probability space (A, (p). Subalgebras {Ai)i£i 
are freely independent iff for any n — 2,3,... and for any Xj £ -^j(j) with £ /, j = 1, . . . , n any 
n-cumulant 

7^„(Xl,...,x„) = 

if there exists a pair fc, I € {1, . . . ,n} such that i{k) ^ 

In sequel we will use the following formula from [TIT which connects cumulants and moments for 
non-commutative random variables 

n k 

(3) ^(Xi...X„) = ^ 7^fc(Xl,X,,,...,X,J [] ^(X,^.+i...X,^.^,_i) 

k=l l<i2<...<ik<n j = l 

with «! = 1 and ik+i = n + 1 (empty products are equal 1). 

The classical notion of conditional expectation has its non-commutative counterpart in the case {A, (p) 
is a W*-probability spaces, that is A is necessarily a von Neumann algebra. Namely, if B C ^ is a von 
Neumann subalgebra of the von Nuemann algebra A, then there exists a faithful normal projection from 
A onto B, denoted by ip{-\B), such that ip(ip{-\B)) = ip. This projection <p(-|i3) is a non-commutative 
conditional expectation given subalgebra B. If X e ^ is self-adjoint then <p(X|S) defines a unique self- 
adjoint element in B. For X e by <y5(-|X) we denote conditional expectation given von Neumann 
subalgebra B generated by X and I. Non-commutative conditional expectation has many properties 
analogous to those of classical conditional expectation. For more details one can consult e.g. [3S]. Here 
we state two of them we need in the sequel. The proofs can be found in [9]. 

Lemma 2.1. Consider a W* -probability space {A,ip). 

• // X G and Y G B, where B is a von Neumann subalgebra of A, then 

(4) ^{■KY)^pMX\B)Y). 

• If \, Z (z A are freely independent then 



(5) 



>p{X\Z) = p{X)I. 
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For any n — 1,2, . . let (xj"'' , • • • , Xq"'') be a family of random variables in a non-commutative proba- 
bility space {Am fn)- The sequence of distributions (/^j^jj(n) ^j) converges as n — > oo to a distribution 
/I if ^J■Q^(.rl) q)^'^^ ^ m(-P) ^ G C (xi, . . . , Xq). If additionally /i is a distribution of a family 

(Xi, . . . , Xg) of random variables in a non-commutative space (y^, <y9) then we say that (X^"-* , . . . , X^"^) 
converges in distribution to (Xi, . . . ,Xq). Moreover, if Xi, . . . ,Xq are freely independent then we say 

that xj"\ . . . , Xg"'' are asymptotically free. 

Now we introduce basic analytical tools used to deal with non-commutative random variables and 
their distributions. 

For a non-commutative random variable X its r-transform is defined as 



(6) rx(z) = £ 7^„+l(X), 



In |36| it is proved that r-transform of a random variable with compact support is analytic in a neigh- 
bourhood of zero. From properties of cumulants it is immediate that for X and Y which are freely 
independent 

(7) rx+Y^rx + rY- 

This relation explicitly (in the sense of r-transform) defines free convolution of X and Y. If X has the 
distribution /i, then often we will write r^ instead rx- 

Another analytical tool is an S'-transform, which works nicely with products of freely independent 
variables. For a noncommutative random variable X its S'-transform, denoted by Sx, is defined through 
the equation 

(8) RxizSxiz)) = z, 
where Rx{z) — zrx{z). For X and Y which are freely independent 

(9) SxY — Sx Sy- 
Cauchy transform of a probability measure ly is defined as 

Jm z-x 

Cauchy transforms and r-transforms are related by 

(10) (^r,{z) + ^^=z. 

Finally we introduce moment generating function Mx of a random variable X by 

oo 

(11) Mx(z)-5] ^(X")z". 

n=l 

Moment generating function and S'-transform of X are related through 



(12) Mx -—Sx{z) 

1 + z 
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3. Free transformations of freely independent free-Poisson and free-binomial 

variables 

A non-commutative random variable X is said to be free-Poisson variable if it has Marchenko-Pastur(or 
free-Poisson) distribution v = i^(A, a) defined by the formula 

(13) V — max{0, I ~ X} 5q + Xv, 

where A > and the measure v, supported on the interval {a{\ — a/A)^, a{l + \/A)^), a > has the 
density (with respect to the Lebesgue measure) 

1 



v{dx) = J'iXa^ - (a; - a{l + A))^ dx. 

2'Kax 

The parameters A and a are called the rate and the jump size, respectively. 

Marchenko-Pastur distribution arises in a natural way as an almost sure weak limit of empirical 
distributions of eigenvalues for random matrices of the form X where X is a matrix with zero mean 
iid entries with finite variance, in particular for Wishart matrices, (see |28j ) and as a marginal distribution 
of a subclass of classical stochastic processes, called quadratic harnesses (see e.g. [T3l). 

It is worth to note that a non-commutative variable with Marchenko-Pastur distribution arises as a 
limit in law (in non-commutative sense) of variables with distributions ((1 — ■^)5q + j^Sa)'^'^ as, N ^ oo, 
see [301 ■ Therefore, such variables are often called free-Poisson. 

It is easy to see that if X is free-Poisson, i^(A,Q!), then 7^„(X) = a" A, n — 1,2,.... Therefore its 
r-transform has the form 

Aa 

^ ' 1 — az 

A non-commutative random variable Y is free-binomial if its distribution (3 ~ I3{a,9) is defined by 
(14) /? = (1 - 't)Io<.<i Jo + /3 + (1 - 0)Io<e<i6i, 

where /? is supported on the interval (a;_, x+), 

x± = 
has the density 

, \ / /,x \/(x — a;_) (x+ — x) , 

Pidx) = ia + 9) ^\ \ L dx. 

27ra;(l — x) 

This distributions appears in CC (unfortunately, the constant a + /3 is missing in the expression given 
for the density part in Cor. 7.2 in CC) as a limit distribution for beta matrices as well as spectral 
distribution for free Jacobi processes in PT| and for a subclass of free quadratic harnesses (see [T^ and 
[11|). The n-th free convolution power of distribution 

pSo + (1 -p)Si/n 

is free-binomial distribution with parameters a — n{l — p) and 9 = np, which justifies the name of the 
distribution (see [34]). 

Its Cauchy transform is of the form (see e.g. the proof of Cor. 7.2 in CC) 

„ , , (g + - 2)2 + 1 - g - ^[{a + 9 - 2)z + 1 - a]^ - 4{1 - a - 9)z(F^ 

(15 Ga^eiz) = — ^ . 

2z(l — z) 

So far the range of parameters a, 9 for which (jl4p is a true probability distribution has not been completely 
described in the literature. In CC the authors seem to assume 9,a > 0, which apparently is not enough 
for a correct definition. On the other hand [17] assumes a,9 > 1 which is too restrictive in general. The 
complete set of parameters is described below. 
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Proposition 3.1. The formula (fT4|) defines correctly a probability (free-binomial) distribution if and only 
if {a, 9) ^ G defined as 

Proof. Recall a result from [llj whicli says that the two parameters Askey- Wilson probability measure, 
which has the form 



(16) .(d^) ^ , (l + a2-2ax)(l + 62-25x)'^^ + ^^^I'^^i'^^ + ^rzr^M>^5._^^ 

is well defined iff a6 < 1. With an additional natural assumption ah ^ 0, this probability law can be 
easily transformed into a free-binomial distribution. Indeed, if we take a random variable X with the 
above Askey- Wilson distribution, then Y defined as 

has a free-binomial distribution (jl4p with parameters 

1 — ah \ — ah 



ip{a,h) 



^a{a — b) ' b{h — a) ^ 
Define now 

(18) H ^ {{a,b) : 1- ab> 0,ah^O}. 

It is rather immediate to see that if (a, h) ^ H then ip{a, h) £ G. 

Conversely, define an equivalence relation on H by (a, b) ^ {—a, ~h) and note that ip is a bijection 
between H/^ and G. 

Finally, referring to the result in [11^ on Askey- Wilson distributions mentioned above, we conclude 
that (ITi|) defines correctly a probability measure (free-binomial distribution) iff (cr, 9) G G. □ 

Note that in "ST the parameters of free-binomial distributions are a ~ n(l — p) 9 = np, and n > 2 
so above conditions are satisfied. Nevertheless, their parametrization does not cover whole G, e.g. the 
situation when one of the parameters is negative, which is allowed. It is worth to note that our derivation 
extends free-binomial distribution to the case a + 9 < 0, in this case continuous part of free-binomial is 
not supported in (0, 1), in case ct < continuous part is supported on {—oo, —1), in case 9 < continuous 
part is supported on (l,oo). 

In CC authors consider complex Wishart matrices N x N, corresponding Gindikin set is {1, 2, . . . ,N — 
2} U [TV - 1, oo) (see [33]). To define Beta matrix as Z = (X -I- Y)-^/'^X{X + Y)-^^^, where X, Y are 
independent Whishart matrices, matrix X + Y must be invertible, so we have pn + Qn > N — 1. For 
existence asymptotic distribution is necessary to pn/N a and qn/N -^• 9 then cr, 6'>0, a + 9 — 1 >0 
so (cr, 9) € G. In this paper further we are concerned only with such case, (cr, 9) G G and cr, 6* > 0. 

The main result of this section, as announced in Introduction, is a direct dual version of Lukacs 
characterization. 

Theorem 3.2. Let (A, ip) be a W* -probability space. Let V and V in A be freely independent, such that 
V is free-Poisson with parameters (A, a) and U is free-binomial with parameters (cr, 9), a -\-9 = \. Define 

(19) X = V^UV^ and Y = V-V5UVi 

Then X and Y are freely independent and their distributions are free-Poisson with parameters {9, a) 
and (a, a), respectively. 

Throughout this paper we use the framework of VF*-probability space, since essentially we work with 
conditional expectations, however the above theorem holds true in a more general setting when {A, ip) is 
a C* -probability space (see [30l Chapter 3]) with exactly the same proof. 
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Proof. Since freeness is defined for subaigebras tfien, witliout any loss of generality, instead of V we can 
take which is free-Poisson with the jump size equal 1. Consider a non-commutative probability space 
{An^ 95 Af) of p-integrable for any 1 < p < 00 random matrices of dimension N x N defined on a classical 
probability space {il, P) with (p{A) — E tr A for any A e An- From the proof of Th. 5.2 and Prop. 
5.1 in CC, in which asymptotic freeness of U and V is stated, it follows that there exist independent 
sequences (Ujv)jv and {'Vn)n of complex N x N matrices such that (Uat, Vjv) converges in distribution 
in the non-commutative sense (as elements of non-commutative probability spaces (Ant^n)) to (U, V) 
which are freely independent. Moreover, XJn is a beta matrix with suitable positive parameters pn, qn 
such that ^ — >■ (T, ^ ^ 9 and Vjv is a Wishart matrix with parameters pn + qN and jjIn, where I^r 
is an TV X identity matrix. It is well known in such case, see e.g. 15 , that 

Xat = v| Vn v| and Yat = Vat - v| Vn Vl 

are independent complex Wishart matrices with parameters {pn, jj^n) and {qn, ^^n), respectively. By 
Th. 5.2 from CC it follows that (X^r, Yjv) are asymptotically free. Therefore, see Prop. 4.6 in CC, 
it follows that there exist freely independent non-commutative variables X' and Y' with free-Poisson 
distributions with jump parameter 1 and rate parameters a and 9, respectively, such that {'K.m,'Vn) 
converges in distribution (in the non-commutative sense) to (X', Y'). 
By asymptotic freeness it follows that 

(20) lim (pAr(P(X„,Yjv)) = ¥'(P(X', ¥')) 

A''— >oo 

for an arbitrary non-commutative polynomial P G C {xi, X2)- On the other hand by the definition of Xjv 
and Yjv 

Vn{PO^n, Yjv)) = ^n{P{^n - V]^/2u^^v]^/^ V]^/'Ua,V]^/')). 
By the tracial property oi lpn the last expression can be written as 

^jv(g(UAr, Va,)), 

for some polynomial Q from C(a;i,a;2). Since (Ujv, V^r) converge in distribution (in non-commutative 
sense) to (U, V) it follows that 

lim (^^(F(X„,Yjv)) = lim (^A,(g(Ujv, Vjv)) = ¥>(Q(U, V)). 
Using the tracial property of we can return from Q to P, so that 

V3(g(u, V)) = if{p{N - v^/^ uvl/^ m^'^)) = ip{P{x, y)). 

Therefore (pH)) implies that for any P e C {xi,X2) 

^(P(X', Y')) = ^(P(X, ¥)). 
Consequently, they have the same distribution. □ 

Corollary 3.3. Let U and V be freely independent random variables in a W* -probability space. Assume 
that V is free-Poisson with parameters 9-\-a and a and U is free-binomial with parameters a and 9. Then 



V2 UV^ = 9{9 + l)a 



2l 
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Proof. By Theorem [O we know that X = UV^ and Y = V - UVi are freely independent. 
Therefore (O of Lemma 12.11 imphes 

(p(X|Y) = ip(X)l 

and 

(^(X^IY) = ip(K^)L 

Due to Theorem 13.21 X is free-Poisson with parameters 9, a. It is well known that for free-Poisson X its 
first two moments are (p(K) ~ 9a and "^(X^) = 9{9 + l)a^. □ 

4. NON-COMMUTATIVE DUAL LUKACS TYPE REGRESSION 

In this section we give the main result of the paper which may be treated as a counterpart of the 
regression characterization of free-Poisson distribution given in Th. 3.2 (ii) of [9 . It is also a non- 
commutative version of Th. 1 of [5] . 

Theorem 4.1. Let (A, if) be a W* -probability space and U, V be non- commutative variables in (A, ip) 
which are freely independent, V has a distribution compactly supported in (0,oo) and distribution ofV is 
supported in [0, 1]. Assume that there exist real constants c\ and C2 such that 

(21) (^(^V-V3UV3 V5UV5^ =ClI 

and 

(22) ^^V- V^UV^)^ V^UV^^ = C2I. 

Then V has free-Poisson distribution, I'iX, a) with A = a + 
and U has free-binomial distribution, l3(a,9). 

/ 1 1 \ " 

Proof. For any positive integer n multiply both sides of (I?!]) and ([2^ by I V2 U V2 j and take expecta- 
tions (fi. Therefore, by Q of Lemma 12. II we obtain, respectively, 

(23) (p(V(VU)") - <y9((VU)"+i) = ci(^((VU)") 
and 

(24) (p(V2(VU)") - 2<^(V(VU)"+i) + (^((VU)"+2) c2(p((VU)"). 
Introduce three sequences of numbers (Q;n)n>o, (/3n)n>o and (7n)n>o as follows 

an^fiiVVD, /3„ = (/.(V(VU)"), and 7„ = (p(V2(VU)"), n = 0,l,.... 



Then equations ([23)1 and ([24]) can be rewritten as 

(25) /3„ - an+l = CiUn 

and 

(26) 7„ - 2/3„+i + a„+2 = C2a„. 

Multiplying (l?5[) and (1^51) by z", z € C, and summing up with respect to n = 0, 1, . . ., we obtain the 
equations 

(27) B{z) - l(A(z) - 1) - ciA(z) 
and 

(28) C{z) - -(B{z) - /3o) + \[A{z) - a^z - 1) = C2^(z), 
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where 

oo oo oo 

A{z) - ^ a„ z", B{z) = J2 Pn 2", C{z) = ^ 7„ z" 

n— n— n— 

and the above series converge at least in some neighbourhood of zero, due to the fact that supports of U 
and V are compact. Note also, that we have A = Myv + 1- 

Before we proceed further with equations ([?7)l and (|28l) we need to establish some useful relations 
between sequences and (7n)- To this end we need to define additional sequence ((5„),i>o, by 

setting 

J„ =^(U(VU)"), n = 0,l,... 
Consider first the sequence (a„). Note that by formula ^ it follows that 

an =7ei(/p(U(UV)"-i) 

+ 7^2[(y5(U)^y3(U(UV)"-2) + (^(UVU)(y5(U(VU)"^-'') + . . . + (y3(U(VU)"-2)^(l[J)] 

+ . . . 

+ 7e„¥'"(u), 

where TZn ~ 7?.„(V) is the nth cumulant of the variable V. Therefore, in terms of (5„'s we obtain 

an = TZl5n-l + 7?.2((5o(5„-2 + 5l5n-3 + • . • + (5n-2(5o) + • • ■ + TIuSq 

and thus for any n = 1, 2, . . . 

n 

a„ = X! X! ■ • • ■ 

Consequently, 

oo n oo n 



Ik 



71 — 1 k—1 ii-\-...-\-ik—n — k n—1 k—1 . . .-\-ik—n—k 

oo oo oo oo 

fe=l ri=fc ii+...+ifc=Ti — i; fc=l m=0 ii + ...+ik=m 

oo / oo \^ 

=1 + ^ 7^fc zM ^ 6,z' . 

k=l \i=0 J 

Therefore 

oo 

A(z) = l + ^7^fc(2A(z))^ 
fc=l 

where A is the generating function of the sequence ((5„), that is A(z) = X^tlo '^i-^^*- Finally, with r being 
the r-transform of V, that is r(z) = J2kLo ^k+iz'^ we obtain 

(29) A(z) = 1 + zA(z)r(zA(z)). 

Similarly, by ([3]), 

/3„ = 7?.ia„ + 7?.2(ao(5„_i + aiSn-2 + . . . + a„_i(5o) + . . . + TZn+iaoSQ 
and thus for any n = 0, 1, . . . 

Tl+l 

A;— 1 2i + . . n — 
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Therefore, 

oo n+l 

n=0 k=l ii + ...+i+fe=n-fe+l 
oo n+1 

= ^ ^ 7?.fc ^ ^ Ctii^ L/Jjx 

n=0 fc=l ii + ...+i+k=n-k+l 

oo oo 

k=l n—k—l ii + .-.+'i+Zc— n — A;+l 

oo oo 

fc=l m— ■ii + .-.+i+fc— m 

oo 

= TZk z^-^ A{z) ^^-\z) = ^(z)r(zA(z)). 
fe=i 

Finally, applying (pUj) we get 

(30) zA(z)r2(zA(z)) + r(^A(^)). 

The formula for 7„ is again based on ([3]) 

+ Tlz[{ao5n-i + . . . + a„_i(5o) + (/3o^o^n-2 + • ■ • + /3„-2<^o)] + • • • + 7?.„+2ao'5o 
and thus it splits in two parts for any n = 0, 1, . . . 

n+2 n+l 

ln = Y '^f' ■■■^ik+Yj'^'' ^^1*^*2 • ■ • ■ 

k—2 i-i+...+ii^—7i — k+2 k—1 ii+.-.+i^— n — /i:+l 

Therefore, also C{z) splits in two parts 

C{z) = Ci{z) + C2{z), 

where 



oc n+2 
ri=0 fe=2 ii + ...+jfc=n-fe+2 



and 

oo n+l 
n=0 k=l ii + ...+ik=n~k+l 

The expression for the second part, C2, can be derived exactly in the same way as it was done for B{z). 
The computation yields 

(31) C2{z) = B(z)r(zA(^)) = zA(z)r3(zA(z)) + r^{z/\). 
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For the first part the derivation is similar though little more complicated 

oo n+2 

oo oo 

= ^7^fcz'^-2^ ^ a,,z''8,,z'^- ...5,, 

k—2 m— i\+...+ik—ra 



Ik 



A{z 

k=2 

Recalling (gH) we get 

(32) = r(zA(z))[r(zA(z)) - U,] + ^(^^(^)) " . 
Finally, ^ and ([321) together with 7^l = /3o give 

(33) C{z) - zA(z)r3(zA(z)) + r\zA{z)) + r(zA(z))[r(zA(z)) - /3o] + ^° 

zA(z) 

Now we can return to the system of equations (|27|) and ([SS)) . Plugging expressions ([29]) and (pO)) into 
([27| we get 

(34) zA(z)r2(zA(z)) + (1 + A(z))r(zA(z)) = ci(l + zA(z)r(zA(z))). 
Define a new function h 

h(z) = zA(z)r(zA(z)) 

and multiply both sides of the above equation by zA(z). Then (I34p it can be written as 

(35) h^{z) = [(1 + ciz)A(z) - l]h{z) + cizA(z). 
Therefore, 

/i3(z) = /i(z) ([(1 + ciz)A(z) - l]/i(z) + cizA(z)) 

and thus 

(36) /i3(z) = ([(1 + ciz)A(z) - if + cizA(z)) /i(z) + [(1 + ciz)A(z) - 1] cizA(z). 
Similarly, plugging ([29l), ^ and ^ into ([28]) we obtain 



zA(z)r3(zA(z)) + 2r2(zA(z)) - /3or(zA(z)) + + 

zA(z) 

2 1 

-- [zA(z)r2(zA(z)) + r(zA(z)) - Po] + ^ [zA(z)r(zA(z)) - aiz] = C2[zA(z)r(zA(z)) + 1]. 

In terms of h the above equation reads 
h^{z) + 2(1 - A(z))/i2(z) - [/3ozA(z) - 1 + 2A(z) - A2(z) + €22"^ A'^ {z)]h{z) ^c^z^A^iz) 

+ /?ozA(z)(l-2A(z)) 

(37) +aizA2(z). 
Inserting (l36l) and ([35]) into (|37|) after cancelations we get 

[(C2 - c?)z + ci + ai - 2/3n]A(z) + /3o - ci 



(38) /i(z) = 



(cf - C2)zA(z) + ci - /3o 
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Plugging it into ([35]) . after canceling A (which is allowed at least in a neighbourhood of zero) we obtain 
the following quadratic equation for A; 

(C2 - c?)[(ci(ai - 2/3o) + C2) + ai + Cl - 2/3o]zA2(z)+ 
+{[(C2 - c?)(ai - ^0) - (ci - /3o)(c2 + ci(ai - 2/3o))]z + (ai + ci - 2/3o)(ai - /3o)}A(z)+ 

+ (/3o-ci)(ai-/3o) = 0. 

The solution of this equation is, 



^ ^'^'^ 2z(c2-c2)(ai+ci-2/3o + (aiCi+C2-2ci/3o)z)' 

where 

Ai = (/3o - ai)(Q;i + Cl - 2I3q), A2 = 2qiCi - aiC2 + C1C2 - ci^occi - Sc^/^o + 2ci^oi 

B = ai + Cl - 2/3o, Co = (ai - 

Cl = 2(ai - /?o) (c2 - c? + ci(/3o - Ci)) , C2 = {c2 " Ci/3o)'. 
We take this branch of square root which makes A holomorphic in 0. It can be rewritten as 

f40) Mz) = 1 + (1 + /^(l - ^) + ^1 (1 - !)) ^ - yi + 2 (ci - + A)) z + (ci + - X)Yz^ 

2/3^(1 + ci(l-i)z) 

where 

c?-C2 C2-C? ci(ai +ci - 2/3o) cl^ ci(q!i + 2ci - 2/3o) „, 

P = -r = 7 A = 2 ^ 2 ^ 2 = 6* + cr. 

ai - Po Cl Cf - C2 C2 - C| Cj - C2 

Note that C2 — c^ is variance so C2 ~ c\ > 0, moreover ci = ip{Y)LpiI — U) > by assumptions about 
distributions, so /? > 0. Similarly, we conclude that a — '^('^)^^(^) > Q. 
Now we use definition of function h, and equations ([55]) . (|40p to get. 



(41) r(zA(z)) = 1 + (^1 - (1 - ^)) ^ - + 2 (ci - /3(1 + A)) z + (ci + P{1 - A))^ 

in a neighbourhood of but not in (note that zA(z) in is 0). On the other hand, r is analytic in 

. , , 1, J fn T., r tn\ r l + (ci-^*(l-A))2-^l+2(ci-/3(l + A))2 + (ci+^;(l-A))2z2 

neighbourhood of 0. Therefore r(0) — hm^-^o 2z — ^ — ~ -^P- 

Consequently, (|41l) uniquely defines r in a neighbourhood of zero. It is easy to check by direct calculation 
that r defined by 

(42) r{z) = \p- ^ 



1 -/3z' 

satisfies (j4T|) and by the analytic extension principle it is unique. Finally we note that (j42|) defines the 
r-transform of the free-Poisson distribution with rate A, and jump size /3. 

It remains to show that V has free-binomial distribution, which can be done through calculating 
S'-transforms. 

Since we know the distribution of U and, as it was noticed before, moment generating function of UV, 
we use equations (|12[) to find corresponding S'-transforms as 

'^u(^) = R\ \ R ^'^^ Sm{z) = — r^- 

pA + pz pA — Cl + pz 

Since U and V are free by © we arrive at 
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Now we use ([U and ([ro]) to find Cauchy transform for V as 
1 



+ f-X+{X-2)z+J(l + f-X+{X-2)zy-A{l-X)z{z-l) 

Gv(z) = ^-7- T ■ 

2z(l — z) 

From (|15p it follows that Gy is the Cauchy transform of frcc-binomial distribution with parameters a, 9, 
(X-f = a,f = e). 

Finally, let us mention that a iy*-probability space with free random variables V and U with, re- 
spectively, free-Poisson and free-binomial distributions, can be constructed in a standard way as a free 
product of two VK* -probability spaces, one containing free-Poisson random variable, second with free- 
binomial distribution. For details see [in],|37]. □ 

Combining Theorems 13.21 and 14.11 we get equivalence, as in the classical situation. 

Corollary 4.2. Let U, V be free random variables in a W* probability space. Then V - Vi/2UVi/2 ana 
Yi/2|[jyi/2 g^j-p jj.pp gyjgf on/y z/ V has free-Poisson distribution and U has free-binomial distribution. 

Proof. The " if part is trivially read out from theorem 13.21 Since freeness implies that conditional 
moments are constant (see (O in Lemma [2TT|) the "only if part follows from theorem 14. II □ 

Acknowledgement. The authors thank M. Bozejko and W. Bryc for helpful comments and discussions. 
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